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1 Introduction 



The understanding of nonequilibrium phenomena owes much to the analysis of simple 
models. For example, kinetic Ising models play an important role in the study of coars- 
ening phenomena [1-3]. The so-called Backgammon model is another example of a simple 
model which brings useful insight into the fields of nonequilibrium statistical physics, slow 
dynamics and aging phenomena [4-6]. It was initially introduced in the context of the 
physics of glasses, as an illustration of the fact that a model possessing entropic barriers 
only could reproduce some of the features of more complex models possessing both energy 
and entropy barriers [7]. For example this model exhibits non-stationary dynamics, with 
aging of two-times quantities [7-11] and violation of the fluctuation-dissipation theorem 
[12]. 

This work is devoted to the behaviour of the Backgammon model at finite temperature. 
We perform an exact analysis of the convergence of the mean energy towards equilibrium, 
and of the relaxation of the fluctuations of the energy in the equilibrium state. We also 
give an analytical treatment of the nonequilibrium dynamics of the model, which becomes 
asymptotically exact at long times and low temperature. This study is adapted from the 
theoretical framework introduced in ref. [11]. It allows us to show that the correlation 
and response functions of the energy are aging at low temperature, to give their scaling 
forms, and to compute the associated fluctuation-dissipation ratio. This analysis is made 
possible by a specific feature of the model, already present at equilibrium, which is the 
existence of a clear-cut separation of slow modes from fast ones, respectively known as ol 
and (i relaxation in the phenomenology of glassy dynamics [6]. A parallel study of the 
dynamics of the density fiuctuations in this model will be the subject of a separate paper 
[13]. 

The model is defined as follows. Consider particles, initially distributed amongst M 
boxes in a known fashion. The energy of a box is equal to if the box is occupied by one 
or more particles, and to —1 if it is empty. The total energy of the system is the sum of 
the energies of each of the boxes; it is thus equal to minus the number of empty boxes. 
There is no reference to a spatial structure in the definition of the model, hence it is of a 
mean- field type. At each time step 1/A^, a particle and a box are chosen independently 
at random. The particle is moved to the box according to the Metropolis rule, i.e., with 
unit rate if the energy does not increase, otherwise with rate e"'^, where /? is the inverse 
temperature. In particular, moves from a box containing one particle to an empty box are 
always allowed, since they do not change the total energy. This dynamics satisfies detailed 
balance. 

The main features of this model are easy to grasp intuitively, and can be explained in 
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simple physical terms. Consider for simplicity the zero-temperature dynamics [7-11]. The 
moves which decrease the energy consist in emptying the boxes which contain only one 
particle, to the benefit of a non-empty box. As time passes, the number of empty boxes is 
ever increasing, and so does the mean occupancy of non-empty boxes. It is therefore more 
and more rare to find a box with only one particle. The dynamics is thus slowed down by 
purely entropic effects, since there are no energy barriers. 

Finally we recall that the present model is similar to a coarsening system which never 
equilibrates [1, 2]. The mean occupancy of the non-empty boxes plays the role of the 
mean domain size in the coarsening system. At zero temperature it grows approximately 
as A{t) ~ Int, while the normalised two-time correlation function of the energy has the 
scaling form [11] 

c(M).©''^;^. (1.1) 



2 Physical quantities and their dynamical evolution 

Consider a finite system, made of M boxes containing particles. Let Ni(t) be the 
occupation number of box number i at time t, i.e., the number of particles contained in 
that box. We have 

M 

J2Ni{t) = N. (2.1) 

i=l 

The model is characterised by the following Hamiltonian 7i, and action S, at inverse 
temperature /3: 

M 

We consider the system in the thermodynamic limit with a fixed density (M, AT — > oo, 
p = N/M fixed). In the following we will be interested in quantities which involve the 
occupation number of one box only. Therefore, taking box number 1 as a generic box, 
we describe the dynamics of the system in terms of Ni{t) only. We denote by fk{t) the 
probability that box number 1 contains k particles at time t: 

fkit) = PToh{Niit) = k}. (2.3) 

Equivalently, fk{t) represents the fraction of boxes containing k particles at time t. 

Restricting the analysis to the homogeneous nonequilibrium initial condition where there 
is one particle per box, we have Ni{t — 0) — 1 for all i, hence p = 1 and 

//c(O) = Sk,i. (2.4) 
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The occupation probabilities fkit) obey the following dynamical equations, whose deriva- 
tion is given in Appendix A: 



dfkit) _ k + 1^^^^^^^ ^ - f 1 + ) fk{t) {k > 2), 



dt A{t) "'"^'^ ' ' V Ht) 

d/i(t) 2 



dt Ait) 
dfoit) 



f2{t) + fi{t)fo{t)-2Mt), (2.5) 



dt 



= flit) -flit) Mt), 



where 



^ = l + (e-^-l)/o(^), (2.6a) 
Mt) = e-^ + (l-e-^)/i(t). (2.6b) 

We introduce the compact notation 

^ = J2^'^dfoit)Jiit)]Mt), (2.7) 

where the matrix A^[/o(t), /i(t)] satisfies 

$^-Mfe4/o(0,/i(0] = 0. (2.8) 

fe>0 

The dynamical equations (2.5) are non-linear, since the transition rates involve A(t) and 
nit), and thus depend on /o(t) and /i(t). We notice that the equation for /o(t) can be 
recast as d/o(t)/dt = /i(t)/A(t) - e-^/o(t). Finally, eqs. (2.5) ensure the conservation of 
the moments 

J2fk{t) = 1, (Niit)) = Y^^fkit) =p = l. (2.9) 

k>0 k>l 

The dynamical equations (2.5) have a simple interpretation. If Niit) = k > 2, the rate 
at which a particle enters box number 1 is unity, while the rate per particle for leaving 
that box is 1/A(t). Hence the total departure rate is fc/A(t). Therefore, if < A(t) the 
occupation number of the box has a tendency to increase, while if A; > A(t) the tendency 
is opposite. The occupation number A^i(t) can also be viewed as the position at time 
t of a random walker on a semi-infinite chain, taking values Niit) = k = 0, 1, . . ., and 
obeying the master equation (2.5). The walk is biased to the left if A; > A(t), to the right 
if A; < A(t), so that the walker is confined around k = A(t) by a restoring force. Hence 
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fkit), the distribution of Ni{t), is expected to be peaked around k = A{t) for long times, 
and for A{t) large enough, i.e., at low temperature. Different rules hold at A; = and 1. 
In particular k = plays the role of an absorbing barrier [9] . 

Note that the average position of the walker {Ni{t)) is fixed by eq. (2.9), so that 
the distribution fk{t) will in fact be bimodal at low temperature. This property will be 
illustrated by the case of the equilibrium distribution (sec Section 5). Finally, there are 
two time scales in this system. On the shorter time scale, A{t) may be considered as 
constant, so that there is equilibration in a fixed confining potential. On the longer scale, 
the potential itself evolves slowly, as A{t) increases. These time scales will be identified in 
Sections 5 and 6 with those of a and (3 relaxation, found in glassy dynamics. 

• Mean energy 

The first quantity to be studied is the energy of the system. According to eq. (2.2), the 
energy of box number i is Ei{t) = — (^iVi(t),o- The mean energy per box of a thermodynamic 
system thus reads 

E{t) = {E,{t)) = -fo{t). (2.10) 
It can be expressed in terms of A{t) as 

1 - Mt) 

(1 - e-^)A(t) 

At low temperature, we have E(t) ^ — 1 + 1/ A{t) up to an exponentially small correction, 
proportional to e~^. The mean occupancy of the non-empty boxes is equal to 1/(1 — 
/o(t)) ^ A(t), again up to an exponentially small correction at low temperature, hence the 
interpretation of A{t) as a characteristic domain size, in analogy with coarsening systems. 

• Energy correlation function 

The correlation function c{t, s) between the values of the energy Ei of box number 1 at 
the two times s (preparation time, or waiting time) and t (observation time) is defined as 

c{t,s) = {E^{t)E^{s)) ~ {E^{t)){E^{s)) , (2.12) 

with < s < t. Note that the present study is restricted to the local correlation function 
c(t, s), disregarding the non-diagonal correlations between two different boxes. 

Since we have {Ei{t)Ei{s)) = Prob{A^i(t) = 0,A^i(s) = O}, we are led to introduce 
the joint probabilities Prob|A'^i(t) = k,Ni{s) = O}. Equivalently, we will consider the 
conditional probabilities gk{t, s) that box number 1 contains k particles at time t, knowing 
that it was empty at the earlier time s: 

gk{t,s) = Prob{iVi(t) = k\Ni{s) = 0}. (2.13) 



m= ^/ (2.11) 
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These quantities obey the dynamical equations (see Appendix A) 

= J2^Mt)Ji{t)Mt,s), (2.14) 

£>0 

with the initial condition 

gk{s,s) = Sk,o- (2.15) 
Eq. (2.14) is now linear in the unknowns gk{t,s). It implies 

J29k{t,s) = l (2.16) 

fe>0 

at all times t > s. 

The energy correlation function then reads 

c{t,s) = fo{s){go{t,s)-fo{t)). (2.17) 

We also introduce the normalised correlation function 

Cit,s) = ^ = '«^l''^-/f\ (2.18) 
c[s,s) l-/o(s) 

such that C(s, s) = 1. 

• Energy response function 

The local response function r{t, s) is a measure of the change in the energy of box number 
1 at time t, induced by an infinitesimal change of the local temperature of the same box 
at the earlier time s. 

More generally, assume that box number 1 is subjected to an arbitrary time-dependent 
inverse temperature Pi{t), while the other boxes (z = 2, . . . , M) are subjected to a uniform 
and constant inverse temperature To leading order at large M, the occupation prob- 
abilities of the boxes i = 2, . . . , M are still given by the fk{t), while box number 1 has a 
different distribution, depending on /?i (t) , denoted by f^^ (t) . 

The dynamical equations for the /^^ (t) are slightly different from eqs. (2.5) (see Ap- 
pendix A): 



d/f^(0 _k + l 
dt A{t) 

dt A{t) 



fkUit) + ftiit) - (l + x|)) fk'it) {k > 2), 

fi^{t)+fi+{t)fS'{t) - {l + ^^_{t))ft{t), (2.19) 



dt 



/._(t)/r(t)-/.+ (t)/o^Ht), 
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with the initial value 



In eqs. (2.19) we have set 



/r(o) = <^M- 



(2.20) 



f^+it) = (1 - flit)) e-^^(*) + hit) Wim - P), 
//_ (t) = 1 - foit) + foit) >V(/3 - piit)), 



(2.21) 



where 

WiAS) = min(l, e"^^) (2.22) 

is the Metropolis acceptance rate associated with a change of action AS (see eq. (A.2)). 
The response function is defined as 



rit,s) = 



mit)) 



Defining more generally 



we have 



hkit,s) 



SPiis) 
5 ft it) 



(2.23) 



/3i(s)=^ 



5Piis) 



f3i(s)=f3 



r(t, s) = hoit, s). 



(2.24) 



(2.25) 



The dynamical equations for /ifc(t, s) are obtained by taking the functional derivative of 
eqs. (2.19) with respect to Pi for /3i(s) = (3. This yields 



^^^^ = Y.MM[foit),f\it)]heit,s), 



(2.26) 



£>0 



with the initial condition 



Therefore 



hkis, s) = (4,0 - 5k,i) /x(s)/o(s). 



J2hkit,s) = 

k>0 



(2.27) 



(2.28) 



at all times t > s. It is to be noted that the above equations are not affected by the fact 
that the derivative of the Metropolis acceptance rate (2.22) is not well-defined at AS' = 0. 
Indeed, d>V(AS ^ 0+)/dAS = -1 and d>V(AS ^ 0-)/dAS = do not coincide. 
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We also introduce the normalised response function 



= 4^ = 44^, (2.29) 

r{s,s) //(s)/o(s) 



such that R{s, s) = 1. 



• Fluctuation-dissipation ratio 

The fluctuation-dissipation ratio X{t,s) provides a measure of the departure of the 
system from equilibrium [6, 14, 12]. In the present case it is defined as 

r{t,s) = X{t,s)^^^. (2.30) 

This definition contains no explicit temperature dependence, in contrast with the fiuctua- 
tion-dissipation ratio of e.g. density fiuctuations. 

The energy fiuctuation-dissipation ratio can be expressed in terms of quantities intro- 
duced above. Indeed dgk{t, s)/ds obeys the same dynamical equations as Qkit^ s) or /ifc(t, s) 
for t > s, since s only enters these equations as a parameter. Moreover, let us integrate 
eq. (2.14) for the gk{t, s) to first order in e = t — s: 

gk{t = s + e,s) = 5k,o + {Sk,i - Sk,o) l^{s)e + O {e^) . (2.31) 

Hence 

^^^ = (4,o-4,i)M«)- (2.32) 
Comparing this last expression with the initial condition (2.27) for the hk{t, s) gives 

hk{t,s)^Us)'^^^^^. (2.33) 



As a consequence. 



= hoit, s) + ^{goit, s) - fo{t)), (2.34) 



so that the fiuctuation-dissipation ratio X{t, s) reads 



or else 



X{t, s) = -^^^^ , (2.35) 

ho{t,s) + ^^{goit,s)-fo{t)) 
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We can check from this expression that X{t, s) = 1 at equihbrium, recovering thus the 
fluctuation-dissipation theorem 

HM)=^. (2.37) 
while < X(t, s) < 1 out of equihbrium. 

3 Generating functions and integral representations 

It is possible to obtain integral representations of the solutions of eqs. (2.5), (2.14), 
(2.26) by means of the generating functions 

F{x, t)^Yl f^^^) ^' ^) = E 9k{t, ^' = E ^) ^"^^ (3-1) 

fe>0 k>0 k>0 

These functions obey the following partial differential equations 

^F(x,t) = (x - 1) (^F{x,t) - ^^F{x,t) - Yf{t)^ , (3.2a) 

|g(x,m) = (^-i) (g(^,M)-^|^ 

d_ 
dt 



G{x, t, s) = ix- 1) [G{x, -Q^G{x, t, s) - Yg{t, s)j , (3.2b) 

H{x, t, s) = {x- 1) (^H{x, t, s) - ^-^H{x, t, s) - Yhit, s)^ , (3.2c) 



where 



Yfit) = (1 - e-^)/o(t) = 1 - 1/A(t), (3.3a) 
Yg{t, s)={l- fi{t))go{t, s) + {l- 1/A(t))^i(t, s) 

= {l-e-f){il-Mt))goit,s) + fo{t)gi{t,s)) (3.3b) 

= (1 - e ^) {go{t, s) + fo{t) — — d^^o*^^' I ' 

Similar expressions for Yh{t,s) are obtained by replacing g hj h in the above three ex- 
pressions for Yg. The initial conditions for eqs. (3.2) are derived from eqs. (2.4), (2.15), 
(2.27): 

F{x,0)=x, G{x,s,s)^l, H{x,s,s) = {l-x)ii{s)fo{s), (3.4) 
while the conservation of moments expressed by eqs. (2.9), (2.16), (2.28) implies 

F(l,t) = ^F(l,t) = l, G{l,t,s) = l, H{l,t,s) = 0. (3.5) 

9 



Eqs. (3.2) can be formally solved by the method of characteristics (see Appendix B), 
yielding 

F{x, t)=(l + {x- l)e-^(*)) e(^-i)^(*'°) + du lC{x, t, u)Yf{u), (3.6a) 
G{x, t, s) = e(^-i)^(*'«) + f du lC{x, t, u)Yg{u, s), (3.6b) 

J s 

H{x,t,s) = iJ,{s)fo{s)IC{x,t,s)+ / dulC{x,t,u)Yh{u, s), (3.6c) 

J s 

with the definitions (B.4), (B.7): 

r{t) = f 
Jo 



du 

D{t,u)= I d^;e^('^)-^(^\ 

J u 



(3.7) 



lC{x,t,u) = (1 _a;)e^H-^W+(«^-i)-f^(*.«). 
Note that, in agreement with eq. (2.33), we have 

u( + \ f / .dG{x,t,s) 

H[x,t,s) = fo{s) — . (3.8) 

Setting a; = in eqs. (3.6), we obtain 

/o(t) = - e-^W) e--^(*'«) + du}C{0,t,u)Yf{u), (3.9a) 
go{t,s) =e-^^^'''> + f dulC{0,t,u)Yg{u,s), (3.9b) 

J s 

ho{t,s)^li{s)fo{s)IC{0,t,s)+ [ dulC{0,t,u)Yh{u,s). (3.9c) 

J s 

Eq. (3.9a) is an implicit non-linear integral equation for A(t), since fo{t), T{t), D{t,u), 
and /C(0,t,tt) are defined in terms of A(t) by eqs. (2.6a), (3.7). Let us finally point out 
the central role played by the kernel 



/C(0, t, U) = er{u)-r{t)-D{t,u) ^ 

ou 



in the following sections. 
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4 Infinite-temperature behaviour 



At infinite temperature, the energy of configurations plays no role in the dynamics. The 
rate at which a particle enters any box is constant and equal to unity, while the rate at 
which a particle leaves the box is equal to unity, per particle. Hence A{t) = fi{t) = 1, and 
eqs. (3.6) provide explicit solutions. 

• Mean energy 

Since Yf{t) = 0, eq. (3.6a) yields 

F{x, t) = {l + {x- l)e-*) exp {{x - 1)(1 - e"*)) . (4.1) 
The occupation probabilities fk{t) thus read 

hit) = ((1 - e-^^ + ke-') ^ -■ (4-2) 

In particular, the mean energy 

E{t) = -fo{t) = (e-* - 1) exp (e"* - l) = -1 + Ig-^* + • • • (4.3) 

converges to its equilibrium value £^eq = — with a relaxation time tiq = 1/2. 

• Energy correlation and response functions 
Similarly, noticing that Yg{t, s) — Yh{t, s) = 0, we get 

G{x, t, s) = exp {{x - 1) (1 - e^"*)) , 

H{x, t, s) = (1 - a;)(l - e"") exp (e"* - 1 + s - t + {x - - e""*)) . 
We thus have 

(l_e^-*)fcexp(e^-*-l) 



9k{t,s) 



k\ 



hk{t, s) = (1 - e"*) exp (e"^ + e^"* -2 + s-t) 



k\ 



(4.4) 



(4.5) 



Note that the distributions (4.2), (4.5) thus found are all related to Poisson distributions. 
In particular the gk{t, s) follow an exact Poisson law, with parameter 1 — e*~*. 

The energy correlation and response functions then read 

c(t, s) = {l- e-«) exp (e-* - 2) (exp(e*-*) - (1 - e"*) exp(e-^)), 
r(t, s) = (1 - e-*) exp(e-" + e""* -2 + s-t), 
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so that the fluctuation-dissipation ratio is given by 

(1 _ e-^)e^"* 

^ (1 - e-^)e^-* + e-2«(l - (1 - e"*) exp(e-* - e^"*)) ' ^^'^^ 

The equihbrium correlation and response functions only depend on the time difference 
9 = t-s: 

ceq(^) = (exp(e-^) - 1) , req(^) = exp(e-^ - ^ - 2), (4.8) 
and they obey the fluctuation-dissipation theorem 

r^m = (4.9) 

The convergence of c{t, s) and r{t, s) to their equilibrium values Ceq(^) and req(^) is in 
e~^^, exhibiting again the relaxation time ti]^ = 1/2, while the latter quantities fall-off 
as e~^, with a relaxation time — 1. The occurrence of these two different relaxation 
times will be explained in Section 5. 

5 Equilibrium properties and convergence to equilibrium 

At finite temperature and for long enough times, the model reaches a stationary state, 
corresponding to thermal equilibrium. In this section, we analyse the equilibrium behaviour 
of the quantities defined above (mean energy, correlation and response functions). In 
particular we compute the relaxation times tiq and tiq-* which characterise the convergence 
of the mean energy to equilibrium, and the relaxation of the local fiuctuations of energy 
at thermal equilibrium, respectively. 

• Mean energy 

The stationary solution of the dynamical equations (2.5) reads 



(/fc)eq = e-^-^ (fc>l). 



(/o)eq = 



Aeq - 1 + e"^«<^ e^"^"'' 



(5.1) 



-^eq -^eq 



where the value Aeq of A(t) at equilibrium is related to temperature by 

e^ = l + (Aeq-l)e^-. (5.2) 

The equilibrium thermodynamics of the model, summarised in Appendix C, shows that 
Aeq is to be identified with the fugacity. 
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The occupation probabilities at equilibrium {fk)eq for A; > 1 are thus Poissonian with 
parameter Aeq, with the exception of {fo)eq, reflecting that empty boxes are energetically 
favoured. At low temperature, we have Aeq ~ ^ 1, so that the {fk)eq form a strongly 
bimodal distribution, with many empty boxes, represented by the weight (/o)eq ~ 1 — 1/Aeq 
at /c = 0, and few occupied boxes, represented by the weight l/Agq concentrated in a narrow 
region of width ^/A^ around k = Aeq. 

• Relaxation time of the mean energy 

The relaxation time 4q of the mean energy can be derived as follows. First, using eq. 
(3.10), we recast eq. (3.9a) into the form 

/o(t) = l-/C(0,t,0)- fdu ^^^ff . (5.3) 

Jo ^^\'^) 

Next, in order to linearise this expression around equilibrium, we set 

A{t) = Aeq + 5A{t). (5.4) 
We have then, to first order in 5A{t), setting s = t — u, 

Mt) = (/o)eq + Sfoit), T{t) = + ST{t), 

Aeq (5.5) 

D{t,u) = Aeq(l -e-^/^«^) +5i:»(t,u), /C(0,t,'u) = K(£, Aeq) + 5/C(0, t, u), 
where 



Aea - 1 + e- 



-A, 



eq 



(5.6) 



Aeq(Aeq - 1) 

SD{t,u) ~ -r^ / dC(e-^/^^<^ - e-'/^^^)6A{t - Q, 

Aeq Jo 

and where the equilibrium kernel K{e, Aeq) = limt_^oo •'^(0, t,t — e) reads 

K(£, Aeq) = exp (-£/Aeq - Aeq(l - C-^/^-)) . (5.7) 

This expression has the following simple relationship to the occupation probabilities {fk)eq- 

K(£,Aeq) = ^MA)eqe-^^/^-. (5.8) 

k>l 

For long enough times, the linearised form of eq. (5.3) reads 

Sfo{t)^^[ duK{t-u,Aeq)6A{u)- ^ [ du6IC{0,t,u). (5.9) 

A-ea Jo Aeq Jq 
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Using eqs. (3.10), (5.6), we obtain 



W) ~ / d£ (i^(£, Aeq) - e-^--=/^- ) SAit - e). (5.10) 

Agq Jo ^ ' 

Taking the Laplace transform of eq. (5.10) leads to the characteristic equation 



=K(p,Ae,)-^— -, (5.11) 



Aeq-1 P+l/Aeq' 

where the Laplace transform K(p, Agq) of i^(£,Aeq) reads 

kip, Aeq) = / d£ e-P^ K(£, Aeq) = Aeq / d^ ^^'^^ e^-^^-^) 

Jo 

^ 1 A'^-i 



(5.12) 



-P+/C/Aeq (/c-1)! 



Notice that the subtracted term in the right-hand side of eq. (5.11) is just the first term 
{k — 1) in the last expansion of i^(p, Aeq). 

Eq. (5.11) has an infinite sequence of real negative solutions, which we denote by 
p = —p^^\ with k > 2. The characteristic values p^^}^ represent the inverse characteristic 
times of the relaxation of energy to its equilibrium value. In particular, the relaxation 
time tiq is the inverse of the smallest characteristic value: 

4q^ = 4)- (5-13) 
P2 

At infinite temperature (Agq = 1), the characteristic values are the poles of K (p, Aeq) 
(except the first one, because of the subtraction), hence p^^^ = k with > 2, and tiq =1/2. 
As temperature decreases, the spectrum of characteristic values is continuously deformed, 
and the solutions of eq. (5.11) stay within a boimdcd distance from the poles of K{p, Agq), 
namely p^j}^ ~ — fc/Agq. At low temperature, the relaxation time ti]^ is exponentially large 
in Aeq. Indeed, we have the expansion around p = 

k{p, Aeq) = 1 - e-^- - Meqe-^-/(Aeq) + ■ ■ ■ , (5.14) 

where [11] 

■'^ n>l e>o 
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By inserting the expansion (5.14) into eq. (5.11), we obtain 



^ ^^^/(Aec), (5.16) 

eq 



I.e., 

n^ e'^^'^ / 1 \ / 21n/3+l \ 

while all the other characteristic times are of order Agq «i at low temperature. 

• Energy correlation and response functions: generalities 

At equilibrium, namely for s ^ tiq"*, the energy correlation and response functions 
become stationary, i.e., invariant under time translation. They only depend on the time 
difference 9 — t — s: 

Ceq(^) = (/0)eq((^0)eq(^) - (/o)eq), (5.18a) 
req(^) = (/lo)eq(^)- (5.18b) 

Their expressions can be derived from eqs. (3.9b, c), which respectively become, setting 
s = t — u in the integrals, 

(^o)eq(^) = e-^-(i-*^"''^'>) + /'deK(e, Aeq)(y,)eq(^ - s), (5.19a) 

(Meq(^) = e-^-K(^,Aeq)+ / ds K (s , A,^) (Yh) ec^iO - s) , (5.19b) 

Jo 

where the equilibrium kernel K(£, Agq) has been defined in eq. (5.7), and where 

{Y,U{9) = (1 - e-^) ((1 - (/i)eq) (^o)eq(^) + (/o)eq(^l)eq(^)) , (5.20) 

and a similar expression for (y/i)eq. First of all, we have the identity 

{ho)eM = -(/0)eq^^^^ (5.21) 

(see eq. (2.33)), which expresses the fluctuation-dissipation theorem (2.37) in the form 

-eq(^) = (5.22) 
generalising eq. (4.9) to any finite temperature. 
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The integral equations (5.19) can be solved by Laplace transforms. Indeed, using the 
definition of and the dynamical equation (2.26) for ho, we obtain 



iho)ec,ip) = (/0)eq(l -P(^0)eq(p)) = . -["^ x ^ • (5-23) 

Aeq - (Aeq - 1) I P + JT^ ^{p, Agq) 

The function in the rightmost side of eq. (5.23) is a meromorphic function, with an infinite 

(2) 

sequence of poles, situated on the real negative axis, which we denote hy p = —pl. , with 

(2) 

k > 1. The pj, give the inverse characteristic times of the energy correlation and response 
functions. In particular the relaxation time of these functions is the inverse of the smallest 
pole: 

4? = ^- (5.24) 
Pl 



req(^) = (/io)eq(^) = J] afee-P^''^ (5.25a) 



We have therefore 

_ (2) 

k>l 

Ceq(^) = (/0)eq((^0)eq(^) - (/o)eq) = E ^e-^^^ (5.25b) 

k>l Pk 

(2) 

where the ak are the residues of the rightmost side of eq. (5.23) at the poles p — —p\. . 
The values of these functions at 6* = read 

'^eq(O) = Y^ak = e-^"'^ , (5.26a) 
fe>i 

\^ «^ <f \ n <f\ ^ (1-e-^-'') (Aeq - 1 + e"^-'') 

Ceq(O) = 2^ -(^ = (/o)eq(l " (/o)eq) = ^ -■ (5.26b) 

k>lPk eq 

At infinite temperature (Aeq = 1), the poles of eq. (5.23) coincide with those of 
i^(p, Aeq), hence p^. = k ioi k > 1, and teq = 1- Again, as temperature decreases, the 
whole spectrum is continuously deformed, and the poles stay close to those of i^(p, Aeq): 
Pj^ PS! — /c/Aeq. At low temperature, the relaxation time teq is exponentially divergent, 
while all the other characteristic times remain of order Aeq. Indeed, by expanding the 
denominator of the rightmost side of eq. (5.23) to first order in p, and using again the 
expansion (5.14), we obtain 

4? - ,2^^'^°^ (AeV-^-/(Aeq) + 1 - Aeq) , (5.27) 
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I.e., 



.n^ 2e^^^ / 2 \ 2e^ / 21n/?+l , 



eq 

(2) 

Let us remark that the relaxation time teq of the energy correlation and response functions 
is roughly twice as large as the relaxation time tiq of the mean energy. Indeed the ratio 
ieqV^eq is equal to two, both at infinite temperature (4q = 1/2, = 1) and at zero 
temperature. We have in the low-temperature regime 

This ratio takes its maximal value 2.20359 for /? = 4.1986. Figure 1 shows the dependence 
on temperature of and 4q\ together with their low-temperature estimates (5.17), 
(5.28). 

We have thus met two spectra of inverse characteristic times, the p^^^ , associated with the 
convergence of the mean energy to its equilibrium value, and the p), , associated with the 
energy correlation and response functions at equilibrium. These spectra have a very similar 
dependence on temperature. Their low-temperature behaviour is characterised as follows: 
the relaxation times ti]^ and become exponentially large in Aeq or in /3, while the 
sub leading characteristic times remain of order Aeq. In analogy with the phenomenology 
of glassy dynamics, the slow modes corresponding to the time scale 6 ~ ti]^ ~ teq"* ~ e^"'^ 
will be referred to as a relaxation, while the fast modes corresponding to the microscopic 
time scale 6* ~ 1 or 6* ~ Aeq will be referred to as (3 relaxation. This separation of the a 
and P relaxation processes in a very clear-cut way is a remarkable feature of the present 
model. 

• Energy correlation function 

Let us now consider the low-temperature behaviour of the equilibrium correlation func- 
tion Ceq(^). Using again the expansion (5.14), we estimate the residue ai as 

Aen - 1 (2-) Q'^"" Q'^"" /-.I 2 , . X 

oi ~ ^79 — Pi ~ A ~ 1-^;^ T^ + -- - - 5.30 

A,\ AeV-^-/(Aeq) + l-Aeq 2 \^ Aeq A,\ ' ^ ' 

This expression shows that the sum of eq. (5.26b) is overwhelmingly dominated by its 
first term (/c = 1), since the relative difference between the full sum and the first term is 
exponentially small. As a consequence, the energy correlation function at equilibrium is 
very close to being a pure decaying exponential, corresponding to a relaxation: 

Ceq(^) «^ e-^<\ (5.31) 
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with exponentially small corrections at low temperature. 



• Energy response function 

In contrast with the simple behaviour (5.31) of the correlation function, the response 
function req(0) exhibits both a and /3 relaxation processes at low temperature. 

— In the regime of /3 relaxation {9 ~ 1), the response function at low temperature is 
obtained by expanding eq. (5.12) for Agq large but p finite: 

A-feA^) = ^-^-^ + .... ,6.32, 



hence 



and 



^eq(^)~ 7T(l + e-'^), (5.34) 



1 

2' 

up to correction terms of order l/Agq. 

(2) 

— In the regime of a relaxation {9 ~ teq ) , the asymptotic fall-off of the response function 
reads 

req(^) ^ aie-^/*e'' , i?eq(^) ^ a^e^'^^e-^/'^l^ , (5.35) 
where ai has been estimated in eq. (5.30). 

— In the crossover regime (1 <^ 6* -C ^eq )) normalised response function i?eq(^) exhibits 
a non-trivial P-to-a plateau value 



flp,-a.e-^=.-^l-^ + ...j«-^l-^ + ...j, (5.36) 

asymptotically equal to 1/2 at low temperature. 

Figure 2 shows log- log plots of (a) the normalised correlation function Ceq{9) and (b) 
the normalised response function i?eq(^), for various values of the inverse temperature (3. 
The different properties of these two functions, discussed above, are clearly apparent. 



6 Nonequilibrium behaviour at low temperature 



The equilibrium analysis of Section 5 shows that, at low temperature, the model exhibits 
slow modes {a relaxation, characterised by the time scales 4q and t^), well separated 
from fast ones (/? relaxation, characterised by the microscopic time scale ~ 1 or Aeq). In 
this section, we show that this separation of modes at low temperature still holds out of 
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equilibrium. This property is at the basis of the analysis presented below. The latter 
consists in an adiabatic elimination of the fast modes (nonequilibrium j3 relaxation), gen- 
eralising the method used in ref. [11]. We give an analytical treatment of the quantities of 
interest (mean energy, correlation and response functions), which becomes asymptotically 
exact for long times {s and t — s 3> 1, irrespectively of the relative positions of these vari- 
ables with respect to the relaxation times tiq"* and tiq"*) and low temperatures. This regime, 
hereafter called the (nonequilibrium) a regime, encompasses in particular the aging regime 
(s and t — s -C tiq , teq"*), and the regime of convergence to equilibrium (s ^ teq , ^eq"*)- 
We begin the analysis by considering the difference 

= AW- 

for 1 but e = t—u ~ 1. Since A(t) is varying slowly in the a regime under consideration, 
we have dA(t)/dt <^ A(t). Hence it is justified to perform a Taylor expansion of A(v) 
around v = t: 

..,-.., = /:a.(^.,-.^^...>^.f^.... ,0.) 

This expansion makes sense when e is small with respect to the characteristic time of 
variation of A{t). Note that performing a Taylor expansion around the lower bound v = u 
would lead to a result equivalent to eq. (6.2). We have kept the first correction term, 
in order to show explicitly that it is of relative order (dA(t)/dt)/A(t) ^ 1. Indeed this 
quantity will turn out to be exponentially small in A(t) throughout the a regime. 

The result (6.2) can be used to derive the following estimates 

^(t,^) ^ A(t)(l-e-^/^(*)), 
/C(0,t,«)ftiexp(^-£/A(t)-A(t)(l-e-^/^W)) =K{e,A{t)), ^^"^^ 

where K{e,A{t)) is obtained by replacing Agq by A{t) in the expressions (5.12). As men- 
tioned above, the approximate results (6.2), (6.3) hold for t 3> 1 but e = t — u ^ 1. It 
will, however, be legitimate to employ them for an arbitrary time difference e, since the 
quantities D{t, u) — A{t) and /C(0, t, u) decay exponentially with e anyhow. 

• Mean energy 

Considering now the integral equation (3.9a) for A(t), we perform a Taylor expansion of 
Yf{u) near u = t. It can be checked a posteriori that it is consistent to keep the first two 
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terms of this expansion. Neglecting furthermore quantities which are exponentially small 
in t, we transform eq. (3.9a) into the differential equation 

fo{t) « e-^W + Jo{t)Yf{t) + Mt)^^ + ■■■, (6.4) 



with 



Jo{t) = I du ]C{0, t,u) 
Jo 

poo 

^ / deK(£,A(t)) = i:(0,A(t)) = 1 -e"^^*), 

Ji{t) = / dw (m - t)/C(0, t, w) 
Jo 

«-rde.Jf(.,AW)=^te)) 
Jo dp 



(6.5) 



= -A(t)e-^^^^/(A(t)). 

p=0 

By inserting the estimates (6.5) into eq. (6.4), using eqs. (2.6a), (3.3a), and consistently 
neglecting terms of relative order e~^(*) or e~^, we obtain 

^^A{k{t\K^\ (6.6) 

with 

^,A,A.,)._L(i_JA^). 

The evolution equation (6.6) for A(t) describes the relaxation of energy throughout the a 
regime. Let us consider the following cases: 

— At zero temperature, and more generally in the aging regime {t <^ tiq ) , we recover the 
result [11]: 

dA 1 , , 

dF - 7(A)' 

hence 



t 



or else 

. / N , , , In In t — 2 , , 

A(t) !^ Int + lnlnt + — \- ■ ■ ■ (6.10) 

It is natural to choose the constant to in order to fulfil the initial condition A(t = 0) = 1, 
yielding to = —0.59962. Anyhow, to only brings corrections to the scaling law (6.10) which 
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are exponentially small in A{t), i.e., comparable to the corrections terms which have been 
neglected during the analysis. 

— At finite (low) temperature, thermal equilibrium is reached for long enough times. The 
convergence of A{t) towards the equilibrium value Aeq is exponential: 

Aeq-A(t)P.a(Aeq)e-*/*-. (6.11) 

The relaxation time tiq reads 

*S ^ -^jmr^i " %^^(^ec), (6.12) 

a^(A,Aeq)/aA|^^^^^ Aeq 

in agreement with the result (5.16) at equilibrium, while the prefactor a(Aeq) is given by 

«(Aec) _ /-^^ ,J 1 



til) In : 



-eq A 

■'^eq 



= / dA , , - , (6.13) 

7o lA(A,Aeq) Aeq- AT ^ 



yielding, after some algebra, a(Aeq) ~ 1 for Aeq large. 

— When time t is comparable to the relaxation time tiq , A{t) changes over from the zero- 
temperature logarithmic growth (6.10) to the exponential convergence (6.11) towards the 
equilibrium value Aeq. The scaling behaviour throughout this crossover can be described, 
to leading order in 1/A, by approximating eq. (6.7) as 

yl(A,Aeq) ^ A(e-^ -e-^^''). (6.14) 

This equation can be integrated, to leading order in 1/A{t) and 1/Aeq. Using the definition 
(6.12), we are left with the explicit interpolation formula 

t Ap 



^eq 



ln(l -e^(*)-^^M . (6.15) 



Figure 3 illustrates the accuracy of the low-temperature analysis in the a regime. Exact 
values for A{t), obtained by numerical integration of the dynamical equations (2.5), are 
found to nicely agree with the prediction of the evolution equation (6.6), for various (large) 
values of the inverse temperature p. 

• Energy correlation and response functions: generalities 

In analogy with the derivation of eq. (6.4), we transform eqs. (3.9b, c) into the following 
differential equations 

goit, s) ^ e-^W + Lo(t, s)Y,it, s) + L,{t, s)^^^^ + • • • , 

(6.16) 

ho{t, s) « Lo{t, s)Yh{t, s) + L^{t, s) ^ + • • • , 
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with 



nt pt 

Lo(t,s)= / duJC{0,t,u), Li{t,s)= / du{u-t)}C{0,t,u). (6.17) 

J s J s 



These integrals differ from Jo{t) and Ji(t) of eq. (6.4) only through the values of their 
lower bounds. They are therefore given by the expressions (6.5), again up to exponentially 
small corrections. 

By inserting the estimates (6.5) into eq. (6.16), using eq. (3.3b) and similar expressions 
for Y/j, and again consistently neglecting terms of relative order e~^ or e~^, we obtain the 
homogeneous differential equations 

1 ac(M)^ 1 



c(t, s) dt r{t, s) dt 

with 



A(A,Aeq)+AMe-^+ ^ 



S(A,Aeq) = , ^ . ^ Aeql^_ ^g_-^g^ 

^ ^' A (A V^/(A) + 1 - A) ^ ' 

Changing the time variables from s and t to A(s) and A{t), we have the alternative form 

1 dc{t,s) 1 dr{t,s) 



c{t,s) dA{t) r{t,s) dA{t) 

with 



-a(A((),Ae,), (6.20) 



These evolution equations describe the nonequilibrium relaxation of the correlation and 
response functions of the energy throughout the a regime at low temperature, as we shall 
see in more detail. 

• Energy correlation function 

The correlation function c{t, s) can be directly estimated from eq. (6.20) with the initial 
value c(s, s) = /o(s)(l — fo{s)), or C{s, s) = 1. Indeed, we know from the analysis of the 
equilibrium properties that this function varies very little in the P regime (t — s ~ 1). We 
thus obtain 

C{t, s) exp ^- y dA' a{A', Aeq) j . (6.22) 
This result can be recast into a multiplicative scaling law: 

^^'''^~$(A(t),Aeq)' ^^-^^^ 
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with 



$(A, Aeq) = exp ^ dA' a(A', Aeq) j • (6.24) 
An interesting interpretation of eq. (6.23) consists in introducing an effective waiting time 

_ $(A(,s),Aeq) _ 1 

d$(A(s),Aeq)/d5 i?(A(s),Aeq)' ^ " ' 

along the hnes of ref. [6]. The correlation function and the effective waiting time exhibit 
several regimes: 

— At zero temperature, and more generally in the aging regime (s <^ ieq"*), we have 



AV^/ A) + l 1/12 ^ , , 

aA,oo) = ^^^ — T — — r~- 1 + x 7 + , 6.26 

^ ^ A (A2e-^/(A) + 1 - A) 2 V A A^ ' ' ^ ' 



hence the approximately square-root behaviour 



$(A(.), oo) ^ CO (A(.)e^(^)) (^1 - J-^ + . . ^ cos'/^ \ns(l + ^^^^^ + •••)• 

(6.27) 

The absolute prefactor cq depends on the initial conditions, according to 



Co = exp ( / dA 



oo 



A, 







a(A,oo)-- ( 1 + - 



(6.28) 



with the natural choice Aq = 1 yielding cq = 6.23367. The effective waiting time reads 

— At finite (low) temperature and for long enough times, the scaling function $ blows up 
exponentially, as 

$(A(s), Aeq) ^ co6(Aeq)e^/*e^' , (6.30) 

with 



«eff = 4? « ^,.\ , = ^'^^'^ (A,V-^-/(Ae J + 1 - Aeq) , (6.31) 

in agreement with eq. (5.27). The prefactor 6(Aeq) can be given a convergent integral 
expression, similar to eq. (6.13). 
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— In the crossover region, when s is comparable to the relaxation times and , the 
scaling function $(A(s), Aeq) changes over from the approximately square- root behaviour 
(6.27) to the exponential behaviour (6.30). The scaling behaviour in this crossover can be 
described, to leading order in 1/A, by approximating eq. (6.21) as 

This equation can be integrated, to leading order in 1/A(s) and 1/Aeq. We thus obtain 

$(A(s),Aeq) « (^e-^(«) -e-^-)"'^', (6.33) 

implying the estimate ln6(Aeq) ~ Aeq/2 for Aeq large. Accordingly, the effective waiting 
time departs from its linear growth (6.29) in s to saturate to the value of the relaxation 
time . This quantity admits the simple expression 

.eff - -j^, (6.34) 

again to leading order in 1/A(s), throughout the crossover. 

The above results allow us to describe the behaviour of the correlation function in the a 
regime, by inserting into the scaling law (6.23) the estimates (6.27), (6.30), (6.33) for each 
of the times s and t. The absolute prefactors cq cancel out, so that the scaling law (6.23) 
does not depend on the initial conditions, in agreement with the fact that only large times 
enter the expression (6.22). In particular, in the aging regime (s and t <S t^), both $ 
functions have the behaviour (6.27), so that we obtain 

C{t,s)^(^y^"^^, (6.35) 



tJ Int 

recovering the result of the zero-temperature analysis [11], recalled in eq. (1.1). In the 
converse situation {s,t ^ teq ), both $ functions have the exponential behaviour (6.30), 
so that the equilibrium correlation function (5.31) is recovered. 

• Energy response function 

The response function r(t, s) is expected to behave differently from the correlation func- 
tion c{t, s), and especially to exhibit an appreciable variation in the P regime {9 = t—s ~ 1), 
in analogy with its equilibrium behaviour, described in Section 5. 

The variation of r{t,s) throughout the nonequilibrium P regime can be determined by 
considering ho{t,s) as a function of $, and by solving eq. (3.9c) accordingly, by means 
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of a Laplace transform, as in Section 5, up to the following differences: take A equal to 
A(s), and not to its equilibrium value Aeq; consider fo{t) and fi{t) as variable quantities, 
given by the dynamical equations (2.5). With these assumptions, closed equations for the 
Laplace transforms ho{p), hi{p) and Yh{p) can be obtained, yielding after some algebra 

r{t,s)^li{s)fo{s)R{A{s),e), (6.36) 

with 

(A(.)-l)(A(.)V^ + ^(A(.),Aeq)) 

Mfois) ^ ^ —-2 -, (6.37) 

A(s) 

and where R{A{s),9) is obtained by replacing Aeq by A(s) in the expression of the nor- 
malised response function -Req(^) = e^'"'''eq(^) at equilibrium, with req(^) being given in 
eq. (5.25a). The result (6.36) holds at low temperature, and for s » 1 and 9 finite. 

(2) 

The response function exhibits a l3-to-a plateau for 1 <S ^ <^ teq , with a plateau value 

rpi(s) « //(s)/o(s)i2pi(A(s)), (6.38) 

with (see eqs. (5.30), (5.36)) 

RpiiAis)) ^ ^ — jrj-, . (6.39) 

^^^^^ A(s)V^^^)/(A(s)) + l-A(s) ^ ^ 

(2) 

The result (6.38) interpolates between the equilibrium value rpi(s) = ai for s ^ teq and 
the zero-temperature value 

'''''^^^ ~ A(s)2/(A(s))(A(s) V^W/(A(s)) + 1 - A(s)) ^^'^^^ 

(2) 

for s -C teq ■ The plateau value (6.38) is the appropriate initial condition to be inserted into 
eq. (6.18), in order to describe the nonequilibrium a relaxation of the response function. 
We thus obtain 

r{t,s)^rpi{s)C{t,s). (6.41) 

The above discussion of the scaling behaviour of the correlation function C{t, s) applies to 
the response function r{t, s) throughout the low-temperature a regime. 

Figure 4 shows log-log plots of the normalised correlation function and response function, 
for various values of the inverse temperature p. The data again nicely agree with the 
predictions (6.23), (6.41) of the low-temperature analysis in the a regime. 
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• Fluctuation- dissipation ratio 

The above results yield the following description of the nonequilibrium behaviour of the 
fluctuation-dissipation ratio X{t, s). In the P regime, we can approximate ^o(^) s) by unity. 
Setting again $ = t — s, and using eq. (6.36), we obtain 

X(t s) ~ A{sffi{s)fo{s)R{A{s),e) 

^ ' ^ ~ A{sfi,{s)fo{s)R{A{s), 9) + A{A{s), Aeq) ' ^ ' ^ 

(2) 

Hence X{t, s) also exhibits a P-to-a plateau for 1 <S ^ <S tgq , with a plateau value 

X ( Hsffiis)fois)R^i{Ais)) 
^'^ ^ ~ A{sfi^{s)fo{s)R^,{A{s)) + A{A{s), AeJ ' ^ ' ^ 

This result interpolates between the equilibrium value Xp\{s) = 1 for s ^ t^ and the 
value 

^^'^'^ A(.)2e-^(^)/(A(l)) + (A(.) - 1)2 ^ ^ " A(^ + ■ ■ ■ ^ ^ ~ (li^ + ■ ■ ■ ^^-^^^ 

at zero temperature, and more generally for s <^ t^ . 

Finally, the evolution equations (6.18), (6.20) imply that the fluctuation-dissipation ratio 
X{t, s) remains equal to Xp\{s) for all subsequent times {9 = t — s 1). In other words, 
the evolution of the fluctuation-dissipation ratio does not couple at all to the slow dynamics 
of the model in the a regime. Figure 5 illustrates this behaviour. 



7 Discussion 



The present work is a continuation of previous studies [9, 11] of the so-called Backgam- 
mon model, introduced in the context of the dynamics of glasses [7]. We have focused our 
attention on the mean energy of the model, its local (or diagonal, see below) correlation 
and response functions at finite temperature, and the associated fiuctuation-dissipation 
ratio. The analysis of equilibrium properties (Section 5) has revealed a remarkable feature 
of the model, namely the possibility of separating in a controlled way the slow modes (o: 
relaxation) from the fast ones (/3 relaxation) . This property is at the basis of the analytical 
treatment of the dynamics (Section 6), which becomes asymptotically exact in the regime 
of nonequilibrium a relaxation: low temperatures and long times (s and t — s ^ 1). This 
approach consists in eliminating the fast degrees of freedom, corresponding to P relaxation. 
It allows a quantitative description of the aging regime (correlation and response decay as 
a power law, the effective waiting time grows linearly), of the convergence to equilibrium 
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(correlation and response decay exponentially with characteristic time tgq of the equilib- 

(2) 

rium situation, the effective waiting time saturates to teq ), and of the crossover behaviour 
between them. The present method is an extension of the approach developed at zero 
temperature [11], interpreted there as an adiabatic approximation. It is worth noticing 
that the physical picture of a and /3 relaxation was absent from this zero-temperature 
analysis. 

The multiplicative scaling law (6.23) of the energy correlation function in the a regime is 
one of the main results of this work. It is consistent with the description of nonequilibrium 
glassy dynamics given e.g. in ref. [6]. Especially, for any three (long enough) times, we have 
C^tijts) = F(C{ti,t2),C{t2,tz)), with F{u,v) = uv. The associated fixed-point equation 
u = F{u,u) only has the trivial solutions u = and u = 1. This simplistic behaviour 
reflects that the correlation function has no non-trivial plateau value. An investigation of 
the density correlation and response functions, and of the associated fluctuation-dissipa- 
tion ratio, will be the subject of a separate paper [13]. The density correlation function 
does exhibit a non-trivial plateau value [12]. 

In this work we concentrated on the diagonal parts of the correlation and response 
functions of energy, which involve the same box at two different times s and t. These 
diagonal contributions are known to be the leading ones in disordered systems, such as 
spin glasses. In the present model, however, non-diagonal correlations, involving say box 
number 1 at time s and box number 2 at time t, are not negligible in general. Let us take 
the example of equal-time equilibrium correlations of the energy. The diagonal correlation 
(E^) - = (/o)eq(l - (/o)eq) peMs a Contribution C^^^g = /32(/o)eq(l - (/o)eq) ~ to 

the low-temperature specific heat per box. On the other hand, the total specific heat at 
low temperature, given in eq. (C.8), reads Ceq ~ 1. We have therefore Ceq -C C^^^^- This 
demonstrates the role of non-diagonal correlations, which screen the diagonal ones almost 
perfectly in this example. 

The present work can also be extended in another direction. In ref. [9] three variants of 
the Backgammon model, with different dynamical rules, were introduced and investigated 
at zero temperature. The model studied in the present work is model A, the other two 
being models B and C. We wish to emphasise that only model A has the property that the 
spectrum of relaxation times at equilibrium can be split in a controlled way into an a and 
a /3 component. We also recall that only model A possesses entropic traps and a very slow 
zero-temperature dynamics (A(t) ~ Int), while A{t) obeys a power-law growth for models 
B and C [9]. These two properties therefore seem to go hand in hand. Finally, we would 
like to mention a generalisation of model B which undergoes a condensation transition at 
finite temperature [15]. 
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Appendix A. Derivation of dynamical equations 

• The case of the fk, Qk} o-iT'd hk 

The rules defining the dynamics of the model are as follows. For a system of M boxes 
and N particles, during every time step dt = 1/N: 

— a particle is drawn at random, which determines a departure box d. The box d is thus 
chosen with probability TTd = Nd/N, since it contains Nd particles; 

— an arrival box a, diff'erent from d, is drawn at random with uniform probability tt^ = 
1/(M-1); 

— the move is accepted according to the Metropolis rule. 

A configuration of the system is given by the occupation numbers of the boxes: C = {A^i, 
A^25 • • • 5 Nm}, with Y2i = ^- We first determine the transition rates appearing in the 
master equation 

Ap(c) ^J2M{C\ C')P{C') - J2 I <^)PiC). (A.l) 

C'y^C C"^C 

Let us consider the rate M.{C"\C) for definiteness. We use the following notation for 
the initial configuration: C = {A^i, . . . , Nd, ■ ■ ■ , Na, ■ ■ ■ , Nm} = {Nd, Na}- Then the final 
configuration is C" = {Nd — I, Na + I}. The Metropolis acceptance rate reads 

W{Nd - 1, iVa + 1 I Nd, Na) = W{AS) = min(l, e"^^), (A.2) 

where the change in action is 

AS = S{Nd -l,Na + l)- S{Nd,Na) = -p{dN„i - Sn^,o). (A.3) 

The possible transitions to be considered are 

C = {Nd>l,Na>0}, AS = 0, W{Nd-l,Na + l\Nd,Na) 

C = {Nd>l,Na^O}, AS^(3, W{Nd-l,Na + l\Nd,Na) 

C = {Nd^l,Na>0}, AS = -p, W{Nd-l,Na + l\Nd,Na) 

C = {Nd = l,Na = 0}, AS = 0, WiNd-l,Na + l\Nd,Na) 

The transition rate M.{C" \ C) is thus given by 

M{C" \C)^M{Nd-l,Na + l\ Nd, Na) = j^7^dT^aW{Nd -l,Na + l\ Nd, Na). (A.5) 

We now consider the reduced configuration characterised by the occupation number of 
box number 1 alone: C = {Ni{t)} = {k}. The transition rates M{C"\C)=M{k ±l\k) 



= 1, 
= 1. 



(A.4) 
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respectively correspond to the case where box number 1 is the arrival box (a = 1: upper 
sign) and the departure box (d = 1: lower sign). They are obtained by summing eq. 
(A. 5) over all the boxes. In the thermodynamic limit (M, N — > oo, with a fixed density 
p = N/M), we getf 



M{k + l\k) = ^^Yl ^d^^W{Nd - 1, iVi + 1 1 Nd, iVi) = p {k> 1), 

d 

M{k-l\k) = ^^Yl ^i^aW{N^ -l,Na + l\N^,Na) = - /o + e-^/o) {k > 2), 

a 

•^(1 1 0) = ^ E ^d7riW{Nd - 1, iVi + 1 1 Nd, iVi) = /i + e-f'ip - /i), 

d 

•^(0 1 1) = ^ E 7ri7raiy(iVi -l,Na + l\N^,Na) = 1. 

(A.6) 

The rates M.{C \ C')=M.{k \ k±l) are also given by these expressions. 

With the notation Mkiifoit), fi{t)] — M.{k\£), and setting p = 1, we obtain the dy- 
namical equations (2.5), (2.7), i.e., 

^^^Y.-^'^^^fo{t),fi{t)]m. (A.7) 

£>0 



The gk{t, s) obey the similar equation (2.14), since the time variable s plays the role of 
a parameter in the dynamics. Finally, the hk{t, s) obey eq. (2.26), as explained in Section 
2. 

• The case of the f^^ 

The only change with respect to the previous case comes from the fact that the local 
inverse temperature of box number 1 is now Pi (t) . The action therefore reads 

M M 

S = pi{t)Ei + pJ2^i = -M)^N„o - /?5](Jiv„o, (A.8) 

hence 

AS = S{Nd -l,Na + l)- SiNd, Na) = Mn^,o - Mn,,!- (A.9) 

f In this limit, it does not matter whether the boxes which are summed upon are assumed 
to be different or not from box number 1. 
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The possible transitions to be considered are 

C^{Nd>l,Na>0}, AS^O, W{Nd-l,Na + l\Nd,Na)^l, 

C = {Nd>l,Na = 0}, AS = pa, W^(iVd-l,iV„ + l|iVd,iV„) =e-/3^ 

C = {Nd = l,Na>0}, AS = -Pd, WiNd-l,Na + l\Nd,Na) = l, 
C = {Nd = l,Na = 0}, AS = /3a-f3d, WiNd-l,Na + l\Nd,Na)=WiAS). 

(A.IO) 

Pursuing the derivation as above, we obtain the dynamical equations (2.19). 
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Appendix B. Solution of eqs. (3.2) by the method of characteristics 

We recall how to solve, using the method of characteristics, a partial differential equation 
of the form 

^G{x, t, s) = {x- 1) (^G{x, t,s)-^^ ^G{x, t, s) - Y{t, s)^ (B.l) 

for t > s, where Y{t, s) is a given function, and where the initial value G{x, s, s) is also 
given. 

The starting point of the method consists in writing the proportionality 

da; dG 
'^'= {x-l)/A{t) = {x-l){G-Yy ^^-'^ 

The first equality in eq. (B.2) yields 

x-l^ye^^^\ (B.3) 

with 



and where the constant y parametrises the characteristic curves. 
The second equality in eq. (B.2) then reads 

^=ye^(^\G-Y), (B.5) 

which can be integrated by varying the constant. This yields, taking in account the initial 
condition, 

G{x,t,s)^e^''-^^^^^''^G(l + {x-l)e^^'^-^^*\s,s^+ duK:{x,t,u)Y{u,s), (B.6) 



with the definitions ^ 

D{t,u) = f dwe^(^)-^(*), 

J u 

}Cix, t,u) = (1 -a;)e^(")-^W+(^-^)'^(*'"). 



(B.7) 
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Appendix C. Equilibrium thermodynamics 



Consider a system made of M boxes containing N particles. As in Appendix A, we 
define a configuration C of the system as given by the occupation numbers of the boxes: 
C = {A^i, A^25 • • • 5 Nm}, with Y^ - Ni = N. The partition function of the system reads 



{iV,>0,^.iV,=iv} 



The combinatorial weight inside the sum represents the number of microstates correspond- 
ing to each configuration C, where a microstate is defined by assigning in which box hes 
each of the particles. The overall in front of the sum takes account of the fact 

that a global relabeling of the identical particles does not aff'cct equilibrium properties. 
This statistics, used e.g. for ideal gases of classical particles, is usually referred to as 
Maxwell-Boltzmann statistics [16]. 

- — — z^i^''~^ for the 

constraint, we obtain 

(C.2) 



'M,N 



with 



= ^-— =e^+e^-l. (C.3) 

fe>0 



We have in particular Zm,n = /N\ at infinite temperature. 

In the thermodynamic limit {M,N oo), the density p = N/M being fixed, the 
method of steepest descent can be applied to the integral (C.2). The saddle-point equation 
reads ds(2;)/d2; = ps{z)/z. The saddle-point value Zc of z, which is by definition the 
thermodynamical fugacity of the model, is thus related to temperature and density through 

e^ = l + (^e/p-l)e"^ (C.4) 

The free energy of the model is an extensive quantity, and its value per box reads 

— /3F = In s{zc) — plnzc = Zc + {1 — p) In Zc — In p. (C.5) 

Physical quantities can be derived along the same lines. For instance, by restricting the 
sum in eq. (C.l) to the configurations {Ni} such that Ei Ni = N and Ni = k, we obtain 
for a finite system 
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and in the thermodynamic hmit 



fe— 1 

(/fc)eq = pexp(-2;c + /55M)^- 

We now restrict ourselves to the case p = 1, considered in the body of the paper. Eqs. 
(5.2) and (C.4) show that the fugacity Zc is identical to the equilibrium value Aeq of A(t), 
while eq. (5.1) for the occupation probabilities coincides with eq. (C.7). 

The low-temperature behaviour of various thermodynamical quantities, such as the 
fugacity Aeq = Zc-, the internal energy £^eq = — (/o)eq, and the specific heat Ceq = 
—0^ dEeq^/dp, is obtained by expanding the above results for /? — > oo: 



ln/3 + 1 
1 ln/9 \n/3-lnP-l 
, 21n/?-l 31n^/?-51n/?-2 



Aeq = /? - In/? + + ^5;^^ + 



^^^ = -1 + 4 + ^ + --- +•••, (C.8) 



The constant equilibrium specific heat Ceq ~ 1 at low temperature is a somewhat patho- 
logical feature of the model [16], that is shared by several classical systems, such as the 
ideal gas and the harmonic oscillator. 
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Figure Captions 



Figure 1: Logarithmic plot of the relaxation times teq (mean energy) and teq (energy- 
correlation and response functions), against inverse temperature j3. Full lines: exact re- 
laxation times (5.13), (5.24). Dashed lines: low-temperature estimates (5.17), (5.28). 

Figure 2: Log-log plot of (a) the normalised correlation function Ceq(^) and (b) the 
normalised response function i?eq(^) at equilibrium, as given by the exact expressions 
(5.25), against the time difference 9 = t — s, for various values of inverse temperature /?, 
indicated on the curves. 

Figure 3: Plot of A{t) against Int, for various values of inverse temperature indicated 
on the curves. Full lines: exact values, obtained by numerical integration of the dynamical 
equations (2.5). Dashed lines: low-temperature estimates (6.6). 

Figure 4: Log-log plot of (a) the normalised correlation function C{t,s) and (b) the 
normalised response function R{t, s), against the time difference t — s, for a waiting time 
s = 100, and various values of inverse temperature indicated on the curves. Full lines: 
exact values, obtained by numerical integration of the dynamical equations (2.14), (2.26). 
Dashed lines: predictions (6.23), (6.41) of the low-temperature analysis in the a regime. 

Figure 5: Plot of the fluctuation-dissipation ratio X{t,s) against the time difference 
t — s, for a waiting time s = 100, and various values of inverse temperature indicated 
on the curves. Full lines: exact values, obtained by numerical integration of the dynamical 
equations (2.14), (2.26). Arrows: low-temperature estimates (6.43) for the plateau value 
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In C(t,s) 




In R(t,s) 




